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Abstract 

Some consequences of the analyticity of the free energy (pressure) of QCD at finite chemical 
■ potential are deduced. These include a method for numerical exploration of the full phase diagram 

[ by a novel use of simulations at imaginary chemical potential to extract Yang-Lee zeroes of the 

grand-canonical partition function. We make use of, and comment on, CPT symmetries, positivity 
of non- linear susceptibilities and the finiteness of screening lengths. We also comment on the 
structure of zeroes expected for the usual picture of the phases of QCD, following a discussion of 

(N 

J> ' the physics of imaginary chemical potential. 
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Motivated by the complex phase structure of QCD at finite chemical potential, /x, pre- 
dicted by effective theories and perturbative analyses there has been a spate of work 
extending lattice computations to finite /i 

mm 

. Most of these methods extrapolate 
data obtained at /i = out to finite /i using, implicitly or explicitly, a Taylor expansion in 
fi of the free energy. As a result, these expansions cannot be continued beyond the near- 
est phase transition to = 0. Much of the interesting phase structure of QCD then lies 
beyond their scope. This paper explores consequences of the analyticity of the free energy 
(guaranteed by time reversal invariance) as a function of complex ^ to suggest a numerical 
technique for exploring the QCD phase diagram. It turns out that analyticity also has many 
other consequences, some of which are touched upon. 

The plan of this paper is the following. First it is shown how analyticity is guaranteed 
by time reversal symmetry. This material is known to practitioners, but serves to introduce 
notation and several concepts I use later. Next I write down the canonical product repre- 
sentation for the partition function and make contact with the Yang-Lee theory of phase 
transitions. A numerical technique is developed to explore the phase diagram using simu- 
lations at imaginary chemical potential, which goes beyond using Taylor series for analytic 
continuation. The need for a different approach is discussed next, with remarks on the phase 
diagram. 

The partition function of QCD, i.e., SU{Nc) gauge theory, with Nj fiavours of fermions, 
each subjected to a common chemical potential, /i, is 

Nf 

Z{T, ^i) = j DUe-^'^'^^ n Det M(mj, /x, T), (1) 



J 

where M is the Dirac operator, Sg the gauge part of the action and the temperature T enters 



the action through boundary conditions on the fields. The free energy, F and the pressure 
are defined by 

PiT, /x) = - ^F(T, ^^) = log Z(T, ^^). (2) 

At any purely real or imaginary value of the chemical potential, the free energy (or the 
pressure) is purely real Is it possible then for the free energy to be an analytic function 
of a general complex fi7 Recall that this requires the real and imaginary parts of F, Fj. and 
Fi respectively, to satisfy the Cauchy-Riemann equations, implying that they are conjugate 
harmonic functions. This rules out Fi = unless F is constant. 
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The answer is that it is entirely possible to have a harmonic function (i.e., a function 
satisfying Laplace's equation) which is zero along the real and imaginary axes and is non- 
trivial elsewhere. Since all even powers of the complex number fi satisfy these requirements, 
it seems that analyticity of F{fi) (at fixed T) requires that it be even in fi. 

Now recall that the transformation fi — /i can be compensated by time-reversal, i.e., 
by interchanging particles and anti-particles. However, this definition is arbitrary; in the 
absence of CP violating terms in the action, thermodynamics remains unchanged by such a 
relabelling. The physics is therefore unchanged by changing the sign of fi. In other words, 
CP symmetry (same as time reversal, by the CPT theorem) implies that F{fi) is an even 
function of fi, i.e., only even powers of fi appear in the Taylor expansion. Since the whole 
chain of logic above can be inverted, this means that CP symmetry is necessary and sufficient 
for the analyticity of F, and its reality along the coordinate axes of fi in the complexplane. 

The simplest use of analyticity is to expand the pressure in a Taylor series ^, O] and 



to determine the coefficients by direct lattice simulations at /i = |a, |S| or by fitting to 
results obtained at imaginary j^. However, this approach is limited, since it only gives 
information on the phase connected to /i = 0. The first phase transition encountered stops 
the extrapolation, and more machinery is required to get around this barrier. This is what 
we develop in this paper. 

At real fi, the path integral measure in eq. ^ is complex Q and importance sampling 
in a Monte Carlo procedure fails although the integral is real. From the argument above, 
it seems that by grouping together configurations related by CP, it should be possible to 
construct a real weight by summing over each such "CP orbit" . Since CP symmetry forms 
a Z2 group, each orbit consists of exactly two configurations, and the sum of the weights in 
;hese two configurations (for real /i) must be real. The simplest example is the Gibbs model 

This model contains Fermions on a 1-dimensional lattice {Nt sites in the temporal direc- 
tion) and f/(l) gauge fields on the links. Since there are no plaquettes, the gauge action is 
trivial and the path integral can be performed easily. The model is a zero (spatial) dimen- 
sional field theory, and, therefore is the quantum mechanics of a single fermion subjected 
to a constant U{1) field. It is possible to choose a gauge where all links have equal value 
U = exp{i6). With fermion mass m, chemical potential fi and defining g^* = z = exp(/i/T), 
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the Fermion matrix is — 



M 



( m qU 
—U^q m qU 
—if^/q m 

V qU 



-W/q\ 



qU 



(3) 



-W/q m ) 



The determinant is- 



Det M = w^' + l/w^' - 2 cosh -Nm' 



(4) 



where w = qU and ma/2 = sinhm'a/2. The lattice spacing is a = 1/TNt. The partition 
function is 

Z{N,z)= J ^^DetM, (5) 

which is real, although the weight is complex. In this simple model a CP transformation 
induces the mapping U -^U'^ . Then, summing over each CP orbit one gets 



^ Det M = 4 cosh(/i/T) cos(A^t^) - 4 cosh(m72T), 

CP 



(6) 



which is real. In the partition function one should compensate the double counting by 
dividing by a factor of two. This is immaterial for expectation values. 

In a more realistic model, the gauge configurations connected by a CP transformation are 
harder to construct. It turns out to be easier to implement the symmetry transformation 
C (which is the same as PT) which is local and maps each link matrix U ^ U'^ . The T 
part of the transformation can be used as before to prove that the free energy is even in /i. 
Moreover, this symmetry transformation preserves the values of 5*^ and Det M (see eq. ^ 
separately. As a result, for the Taylor expansion of DetM summed over each PT orbit — 



Det M(m, T, /i) = Det M(m, T, 0) 



1 +ai/i+ y/i + -^/^ + 



(7) 



one has 3? = ^3 Oj = for all odd i. In other words, summed over PT orbits, the measure 
is even in /i. 

Since this summation over a symmetry orbit leads to a weight that is the sum of two 
determinants, it is clear that the measure cannot be rewritten using local pseudo-fermion 
fields. Disappointingly, these symmetries turn out to be unexploitable in this form for 
molecular dynamics algorithms paralleling those which are used at zero chemical potential 



Z(T, /i) = Tr exp 



(8) 



IjJ. However, by ensuring analyticity, they lead to other numerical approaches to the 
problem, as I show next. 

From the Hamiltonian expression for the partition function, 

H fiN 

(where H is the Hamiltonian and the number operator) it is clear that Z has a periodicity 
of 2itT in the imaginary part of the chemical potential — 

Z(T, /i) = Z(T, /i + 27iikT) for any integer k. (9) 

Functions with such periodicity are best analysed in terms of the fugacity, z = exp(/i/T). 
Entire functions in /x are regular in the complex z plane without the origin, i.e., the punctured 



z plane 



12| . Then the partition function can be written in the canonical product form 



Z{T,z) 



-F{T,z)/T 



n 1 



ZniT) 



Q 

(10) 



where Zn{T) are the zeroes of Z, and F is the regular part of the free energy. CP symmetry 
constrains the set of zeroes, {zn} to be symmetric under inversions in the unit circle. Branch 
points in the full free energy, F = —TlogZ, then do not appear in the regular part F. The 
periodicity of Z in fi also implies that 

F{T, /i) = F(T, fi + 2mkT) for any integer k. (11) 

The appropriate setting for further analysis is the Yang-Lee theory of phase transitions Q| . 

According to this theory, phase transitions occur at points in the parameter space where 
dense sets of zeroes develop into pinch points. In the present understanding of the phase 
diagram of QCD, one expects pinch points on the real /i axis at small T where interesting 
new phases of QCD appear We shall come back to this after a digression on ideal gases 
which throws more light on the question of periodicity in the imaginary part of /i. 

An ideal gas has the free energy — 

1 



F{T,^i) = N,NfV 



(12) 



180 '6"" ■ 127r2' 
which does not have the required periodicity; no polynomial in /i can. By examining the 
Hamiltonian formulation of an ideal gas, we can show where the periodicity is hidden. The 
free energy can be written as 



F(/x,r) 



N^NfVT^ 
37^ 



+ 



1 



ze^ + 1 / z + 1 



(13) 
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Since z is unchanged under ^ ^ fi + 2niT, this has exphcitly the required periodicity. 
By the substitution y = exp x, the integral can be written in terms of fourth order Euler 



Polylogarithm functions 15 1 



7r2 



Li4 (--) + Lu{-z) 



which also explicitly retains the periodicity. In fact, each of the terms has this periodicity 
explicitly. This property is lost in a Taylor expansion of the polylogs in fi/T — 

_ Ttt^ SCs tt^ 2 log 2 3 x'^ x^ x'^ x^ 
Li4 (-e j - -— -—^-^^ - " 48 " 480 ^ 40320 " 1451520 ' ' ' ^ 

where the remaining terms are odd in the expansion parameter. The reason is that in taking 

the logarithm of the argument one makes a choice of the Riemann sheet, which is precisely 

equivalent to working within a single strip of width 27iiT. When the sum in eq. (I14|l is 

expanded in a Taylor series in /i, it is clear that every term cancels, apart from those which 

appear in eq. p2|) . Periodicity can be restored by averaging this expression over all Riemann 

sheets of the logarithm, the result of which is expressed concisely in eq. ()14p [iG^ . In general, 

interactions modify the integral expression in eq. (fT^ . without spoiling the periodicity. This 

upsets the conspiracy which cancels all higher terms in F, leading to non-vanishing values 

for general Taylor coefficients, i.e., the non-linear number susceptibilities 

We return now to considerations of phase transitions in QCD through the Yang-Lee 

mechanism. The set of zeroes, Zn{T), of the partition function contains the full information 

needed to identify all the phase transitions in the theory through the Yang-Lee mechanism 

jl^ . But for all this to make sense, one needs to check that a thermodynamic limit exists, and 

that, in this limit, the pressure is a continuous and monotonically increasing function of z. 

This is Theorem I of |lJ], and is proven there in the non-relativistic limit when interparticle 

potentials obey two conditions — first that there should be a limit to the number of particles 

that can be accommodated in a finite box, and second that the interactions between them 

)e of finite range. Non-intuitive phenomena can occur when these conditions are violated 



17[. All work on QCD assumes the validity this theorem. However, it is useful to first 
examine the basis for this assumption. 

In the low temperature phase, where the carriers of baryon number are the baryons them- 
selves, the problem is essentially non-relativistic, since the baryon mass is much greater than 
the temperature. There is ample experimental evidence that the inter-nucleon potential sat- 
isfies both conditions necessary for the theorem to hold — the baryon-baryon interaction has 



a hard core repulsion and a short ranged Yukawa interaction at long distances. In the high- 
temperature phase the problem is relativistic, since one expects that the carriers of baryon 
number are the quarks, whose masses are much less than the temperature. Quarks and 
antiquarks are freely created and destroyed, and, as is well-known, the correct reformulation 
of the problem is to work in ensembles where the excess (or deficit) of quarks over antiquarks 
is fixed (the canonical ensemble) or is conjugate to the chemical potential (the grand canon- 
ical ensemble). The second requirement, of finite range of interactions, is straightforward. 
The gauge interaction, which could be long-ranged, is known to be screened in the electric 
sector [l^. In the magnetic sector there is no evidence for screening, but there is evidence 
of confinement and hence a mass- gap 19j. All other screening lengths are smaller in the 
continuum limit ^8]. While these phenomena are best established for = 2 and 3, they are 
expected to hold also for other Nc- 

If we accept these arguments, or otherwise directly assume the first theorem of Q| — that 
the pressure is a continuous increasing function for z > 1, then an immediate consequence 
is that the Taylor expansion coefficients of P{T, z) in z at fixed T are non-negative. Now, 
the generalised susceptibilities [6], are Taylor coefficients of P{T,fi) in an expansion in fi. 
Since the expansion in fi can be obtained from that in z = exp(/i/T) by expanding out 
the exponential, which is itself a convex function, it follows that the sum of all generalised 
susceptibilities of a given order are strictly positive. 

This theorem of Yang and Lee is expected to hold for a continuum theory in the ther- 
modynamic limit of infinite volume. It is also expected to hold for a continuum theory in 
a finite volume much bigger than the range of interactions, even when the limit is taken 
with a sequence of cutoffs, holding the physical volume fixed. On the other hand, at any 
intermediate step in this procedure, the convexity argument may fail if some of the screening 
masses are far from their continuum limits and not sufficiently small. Negative values of 
some of the susceptibilities may then be observed at finite volume. However, convexity is 
regained either by increasing the spatial volume or on taking the continuum limit [2^. 

Further practical use of this theory necessitates estimation of the set of zeroes, {zn(T)}. 
Further information hinges on the connection between the grand-canonical and the canonical 
partition functions. The partition function in eq. (0) generates expectation values in a grand 
canonical ensemble with respect to quark number. A partition function for fixed quark 
number, Q, i.e., in a canonical ensemble has also been used before 

mm. 

One 
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defines this canonical partition function by tlie Fourier transform 

Q(T, Q) = ^ r d(PZ{T, 20) cos Q0, (16) 

where Z{T,i(f)) denotes the grand canonical partition function in an imaginary chemical 
potential fi = icf). One can invert this formula to write 

VNsNjN^ 

Z{T,z)= Y: QiT,Q) ^« + ^ , (17) 
Q=l ^ 

where Ns is the number of components of the Dirac spinor (A^^ = 4 in four dimensions). 
Note that on any lattice with V sites on each spatial slice, Z is a polynomial in z of degree 
VNgNfNc. The Z^^ symmetry of the gauge action can be used to show that Qn = unless 
n is a multiple of A*",, (23^. The remaining coefficients are real and non-negative. Hence 
the VNgNf distinct zeroes of Z cannot lie on the positive real axis of z. The Yang-Lee 
picture now follows, with phase transitions occurring when ^ ^ 00 and the zeroes pinch 
the positive real z axis. Note that the finite Laurent series on the right of eq. (fT7j) is typical 
of a relativistic theory with its CP symmetry. In non-relativistic theories this is just a finite 
polynomial. 

This gives a method for determining the transition line in the (/i, T) plane. Construct the 
partition function at imaginary chemical potential and, by Fourier transforming it, build up 
the coefficients of the polynomial in z. Transform this information into that of the zeroes 
of the polynomial and check whether they begin to pinch the positive real z axis. Since 
this information finds all the pinch points, this is one way to go beyond the present day 
lattice computations and look for signals of phase transitions to the conjectured colour 
superconducting phases. Some results will be presented elsewhere. 

This is an appropriate place to discuss a couple of interesting points about this algorithm. 
An absolute normalisation for Z is not needed in any practical implementation of this 
method. An overall multiplicative factor in Z appears as such in the Q's, and hence is 
irrelevant to the zeroes of eq. (fTTj) . Simulations will never give an absolute normalisation for 
Z . However, combinations of simulations at several different values of can be combined to 
give the relative values of Z, as is done in setting up a multi-histogram reweighting. There 



have been preliminary investigations of methods for constructing Q's 

A second practical question is that of statistical errors. Any simulation will necessarily 
contain statistical errors, and the question naturally arises whether they allow a determina- 
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tion of the pinch point of the Lee- Yang zeroes with any degree of accuracy. Since the pinch 
points he on the real axis, after appropriate scahng, the zeroes in the vicinity of this point 
are the solution of + 1 = 0, with being a large integer of the order of V. The zeroes of 
this polynomial lie at the points z„ = exp[2TTi{n + l/2)/N]. The effect of simulation errors 
is to change this equation to a general polynomial of order k — 

N 

Y^aiz' = 0, (18) 

1=0 

where the coefficients are real and we can write |aAr — 1| = tN, |ao — 1| = and \ai\ = 
(for < z < A^). We model the errors as being independent random numbers drawn from 
Gaussians of width a. For sufficiently large statistics a is small, and the deviations of the 
roots of eq. (fTHj) from Zn, Sn, can be treated as linear in the e^. It is then a straightforward 
exercise to show that = and the RMS error, ^ (5^) — = aZn/N. Thus the 

problem is statistically well conditioned. The estimates of zeroes are unbiased and the 
errors fall inversely with the square root of the statistics and inversely with the system 
volume. 

We turn to the question of phase transitions at purely imaginary yU, since this yields 
further insight into the locations of the zeroes of the partition function. We note at the 
outset that crucial convexity theorems (such as the second law) fail at imaginary /i, and 
hence thermodynamics in its usual sense does not apply. The idea of is to use the 
peak in a response function, x, (i.e., second derivative of the free energy with respect to an 
intensive parameter) to identify putative critical points at imaginary /i and finite volume, 
analytically continue them to real /i, and then take the infinite volume limit. Note several 
subtleties in the argument. First, for a function xif^'^)^ the extrema for imaginary chemical 
potential, /x^ < 0, are obtained by solving for the zeroes of the first derivative. This does 
not necessarily give the extrema at real chemical potential, since the solution of x'if^l) = 
does not guarantee that x'{~f^l) also vanishes. A second, related subtlety can be seen in a 
double Taylor expansion for any response function at finite volume — 

X(T/T„;x/T) = (1 - T/TS , (19) 



where the Taylor coefficients Cnm are independent of T and fx [3". At the critical end-point, 
Te^^J's) this series must sum to the divergent quantity |1 — T/TE{fiE)\~'^ (with 7 > 0) in the 
infinite volume limit. There may also be a critical point, {T'^,iJ,'^) at an imaginary chemical 



potential, (/i^;)^ < 0. Clearly the sum in eq. p9|) differs at the two points, i.e., even if 
x(TE/Tc,fiE/TE) and x{Te/Tc, ^-'e/T^) are both divergent, neither x{TE/Tc,ifJ.E/TE) nor 
xiT^Tc, n'E/iT^) need to diverge 

2 



o 



A 
C 

V 




-2% -37C/2 -7C 



37C/2 27C 



-i|aT/N, 



FIG. 1: Number densities, (n) in the model of 



for the values of Q(T, 0) marked in the figure. 



In the context of our earlier arguments, it is more natural to look for accumulation points 
of zeroes of the partition function for imaginary chemical potential. This leads to the third 
subtlety, and a genuine difference between real and imaginary chemical potential due to the 
lack of convexity mentioned earlier. The theorem of Yang and Lee which we have earlier 
quoted allowed us to express the partition function as a finite (at all finite volumes) Laurent 
series in z with non- negative coefficients feq. ITTj). This constraint on the coefficients ensures 
convexity and prevents zeroes at real positive z, i.e., for real /i. Only in the limit of infinite 
volume can the zeroes pinch the real line and cause a phase transition. However, this same 
structure allows, (and, in some models, may force) the occurrence of zeroes on the unit circle 
in z (i.e., for imaginary /i) even at finite volume. 

Examples abound. We point out the extension of the Gibbs model to SU{Nc) gauge 
group where the partition function has been computed 2^, and found to be 



Z = 2cosh[N^n/T] + smh[{N^ + l)m7T]/ sinhfm'/T]. 



(20) 
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This is of the form in eq. p7|) . with Q{T,Nc) = 2, Q{T,0) being the second term in 
the expression on the right, and all other coefficients vanishing. In the limit m ^ 0, 
Z = 2 cosh[Ncfi/T], and the zeroes of Z are at 



J{2n+l)n/2Nc 



(0 < n < N^). 



(21) 



Note that there are no positive real Zn- The zeroes lie exclusively on the imaginary axis in 
the /i-plane. The nearest zeroes to the positive real z axis are ai z = explain /2Nc). In this 
0+1 dimensional problem there is no phase transition for any finite A*",,. However, a pinch 
point develops at z = 1 as iVc — > cx), very beautifully illustrating van Hove's theorem. 




G 



Re ILL 



Im \i 



FIG. 2: Phase diagrams along the real and imaginary axes in fi for Nc = 3 and Nf = 2 for finite 
quark mass. The two insets show the lines of Yang-Lee zeroes pinching the phase transition points 
on the real z axis. The phase diagram on the left constrains only the pinch points for the lines 
marked a and b, not their shapes and orientations. CP symmetry causes the whole phase diagram 
to be replicated for ^fj, < 0, and the corresponding Lee- Yang zeroes to be replicated inside the 
unit circle (dotted line) by inversion. 

The zeroes on the imaginary fi axis give rise to logarithmic branch point singularities of 
the free energy evaluated at imaginary chemical potential. The quark number density, all 
quark number susceptibilities, the energy density and all its derivatives diverge at all these 
points for all values of A^^^. These divergences are thermodynamically meaningless since they 
occur without the necessity of taking a thermodynamic or continuum limit. At finite quark 
mass, the zeroes remain on the unit circle but move about as long as Q(T, 0) < 1. When 
the mass becomes larger, these zeroes leave the imaginary axis, but they are still visible as 
maxima of these observables (see Figure^, and can be interpreted as crossover points. The 



11 



effect of mass seems to be one of the main differences between the perturbative and strong 
couphng computations of 

These crossovers are generic at imaginary chemical potential — Potts models (or 0{N) 
models) with an ordering field in one direction and an imaginary transverse field show such 
behaviour, as do massive fermions subjected to imaginary chemical potential. Removing the 
first ordering field converts them into branch point singularities. Both pose dangers to nu- 
merical simulations, since these non-thermodynamical singularities are easily misinterpreted 
as strong first order transitions, or developing critical points. If one indeed misidentifies 
such points, then finite size scaling studies are also misleading, since one tries to apply 
inappropriate homogeneity relations arising from thermodynamic considerations. The idea 
of locating the zeroes of the partition function using data collected at imaginary chemical 
potential is to bypass these problems altogether. 

The connection of the results of [3] and the phase diagram of QCD at real /i seems to 
be rather subtle. Here we make a hypothesis which is consistent with the observed absence 

n 

of a critical point with Z3 symmetry j2a] and with present understanding of the physics of 
the Zn^ phases [2^. Figure El is a cartoon of the conjectured phases of QCD {Nc = 3 and 
Nf = 2) with finite quark mass Q]. There are three lines of first order phase transitions 
(labelled A, B and C) in the figure. These meet at a triple point, where the normal phase, 

the plasma phase and a colour-superconducting phase coexist. The line A, separating the 

I 

normal and plasma phases ends in a critical point 28] . The line B is expected to continue to 
infinity. For imaginary chemical potential, only the locus of zeroes of the partition function, 
the lines F and G, are shown. 

There are also two insets showing the positions of Yang-Lee zeroes which might lead to 
such phase diagrams. In this figure only the pinch points are constrained — the remainder 
of the lines 'a' and 'b' are allowed to twist, merge or intersect without changing the phase 
diagram Q. The important potat is that in a range of temperatures, there are two pineh 
points on the real axis [20|. As the temperature rises, the pinch point on the right, labelled 
'b', moves to higher and higher fir, whereas the one on the left, labelled 'a', lifts away 
from the real axis at the critical end point, Te- Its effect can still be felt at the crossover 
temperature, T^, seen in finite temperature simulations at /i = 0. Above the line 'b' 
continues to move outwards, but there remain isolated Nc zeroes on or near the unit circle 
(drawn with a dotted line), depending on the quark mass, which give rise to the crossovers 



12 



(F and G) observable at imaginary fi. 

We end by summarizing the contents of this paper. The complex analytic structure 
of the free energy (pressure) of QCD was investigated. It was shown that CP symmetry 
makes the free energy an analytic function of the chemical potential fi at all temperatures 
T. This symmetry is unexploitable to give an algorithm for a Monte Carlo procedure that 
works through a molecular dynamics algorithm. Exploiting analyticity, it is possible to 
write a canonical product representation of the partition function, and made contact with 
the Yang-Lee theory. It was argued that analyticity, monotonicity and continuity of the 
pressure implies that non-linear susceptibilities of arbitrary order must exist and be non- 
negative at all T. It was further argued that this is connected with electric screening and 
magnetic confinement in high temperature QCD. I showed how to use canonical partition 
functions, derived from Monte Carlo simulations at purely imaginary chemical potential, 
to obtain information on the phase structure of QCD. Finally, a plausible phase diagram 
of QCD was used to anticipate the behaviour of the zeroes that such numerical techniques 
might give rise to. 
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